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1 Introduction 



The Anderson localization is one of the most spectacular disorder-induced phenomena pp. 
The one- dimensional situation is especially well-understood Consider for definiteness 
the Schrodinger equation for an electron moving on a line, in a disordered potential V{x). 
Even in the presence of an infinitesimal amount of disorder, all eigenstates are exponen- 
tially localized, with a localization length ^ = I/7, where 7 is the Lyapunov exponent. 
As a consequence, the typical conductance of a disordered sample falls off exponentially 
with its length L. More precisely, the zero-temperature conductance g of a one-channel 
sample is related to the transmission T across the sample by the two-probe Landauer 
formula [H]: 



The theory of one-dimensional localization predicts that the transmission T is a widely 
fluctuating quantity in the insulating regime (7L S> 1), so that the meaningful quantity 
to consider is its typical (i.e., most probable) value exp((lnT)). The mean (InT) grows 
linearly with the sample length. 



and the ratio (In T)/L is self-averaging, in the strong sense that all the cumulants of 
(InT) grow linearly with T In other terms, the statistics of T (resp. of InT) is 
similar to that of the partition function (resp. of the total free energy) of a disordered 
thermodynamical system. This deep analogy appears clearly in the framework of the 
transfer-matrix formalism, especially for discrete (tight-binding) models jH E] . 

So far, it was implicitly assumed that the energy E of the incoming electron is above 
the mean of the disordered potential (usually taken to be zero). Much less is known 
on the converse situation of tunneling through a disordered barrier, where the mean 
potential {V{x)) inside the sample is non-zero, and higher than the energy E. More 
generally, a tunneling situation is met if the disordered sample is periodic on average, 
and if the energy is in a gap of the underlying average structure 0. 

It has been put forward, seemingly in |6 for the first time, that a weak disorder 
enhances the transmission in such a tunneling situation. In their subsequent work [7j, 
the authors show that a weak disorder increases both the mean conductance (propor- 
tional to (T)) and the mean resistance (proportional to (1/T)) of a tunnel barrier. This 
disorder- induced enhanced tunneling effect is paradoxical, because random impurity po- 
tentials usually lead to additional scattering, which hinders transport. In any case, a 
strong enough disordered potential is expected to have the usual effect of reducing the 
transmission. Putting together these observations, it can therefore be anticipated that 
the transmission reaches a maximum in an intermediate crossover regime, corresponding 
to a moderate amount of disorder. This non-monotonic behavior of the transmission as a 
function of the disorder strength seems to have been overlooked so far (see, however, jH] 
for attempts earlier to E], and [0] for a recent discussion of the effects of a weak 
disorder on gap states). 





(InT) ^ -27L, 
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Our aim is to provide a quantitative analysis of the non-monotonic behavior of the 
disorder-induced enhanced tunnehng transmission. We restrict this study to the regime 
of most physical interest: 

a (correlation range of potential) <^ 1/K (penetration length) <^ L (barrier length). 

(1.3) 

The first inequality implies that the fluctuations of the disordered potential are short- 
ranged, so that the latter can be modeled as a Gaussian white noise. The second in- 
equality implies that the transmission of the barrier is exponentially small, even in the 
absence of disorder. This suggests that (InT) will be the right quantity to consider. 

For completeness, we first give in Section 2 an overview of tunneling through a clean 
barrier. Section 3 then deals with tunneling through a square disordered barrier (the 
mean and the width of the potential are constant), while Section 4 is devoted to the gen- 
eral case (the mean and/or the width of the potential vary smoothly across the barrier). 
A summary and an outlook are presented in Section 5. 



2 Tunneling through a clean barrier 

We begin with a reminder on the well-known problem in Quantum Mechanics ^U] of 
tunneling through a clean barrier. In reduced units [h = 2m = 1), the one- dimensional 
Schrodinger equation reads 

-i)"{x) + V{x)i){x) = Ei){x). (2.1) 



Square clean barrier 

Consider first the simple case of a square barrier of length L. The potential is constant 
in the barrier: 

V{x) =Vq (0 < X < L), (2.2) 

and vanishes elsewhere. 

In the situation of interest, the energy E of the incoming particle is in the range 
Q < E < Vq. The wavevector p of the particle and its inverse penetration length K in 
the barrier read 



p=^/E, K= ^V^ - E. (2.3) 

The reflection and transmission amplitudes r and t are determined by looking for a 
solution to ()2.1|1 of the form 



%Ij{x) = < 



ae^^ + 6e--^^ (0 < x < L), (2.4) 
te'P(^--^) {x>L). 
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Expressing the continuity of iplx) and of its derivative at x = and x = L provides four 
linear equations, whose solution yields 



{p"^ + K^)smhKL 



(p2 - K^) sinh KL + 2ipK cosh KL 
2ipK 

{p^ - K^) sinh KL + 2ipK cosh KL ' 



(2.5) 
(2.6) 



Throughout the following, we will be mostly interested in the transmission intensity 
coefficient (or transmission for short), 

T = |t|2, (2.7) 
which enters the Landauer formula In the present case, (|2.6|) yields 

In the regime ()1.3|) . where the barrier length is much larger than the penetration length, 
the transmission falls off exponentially, as 

r-^-^^^exp(-2irL). (2.9) 

All subsequent results for the transmission will be given with exponential accuracy, in 
analogy with (|1.2|) . Neglecting the pref actor, we thus rewrite (|2.9|) as 



\nT^-2KL. (2.10) 



Arbitrary clean barrier 



Let us now consider tunneling through an arbitrary clean barrier. The potential V{x) is 
larger than the energy E for < x < L, so that the inverse penetration length reads 



K{x) = ^V{x) - E. (2.11) 

We assume that the potential has a smooth profile across the barrier, i.e., the length 
scale over which V{x) or K{x) varies is of the order of the barrier length L itself. 

The transmission can be determined along the lines of the previous case, by seeking a 
solution to 1)2.11) of the form ()2.4|) . where exp{±Kx) are replaced by the two elementary 
solutions u{x) and v{x), with initial values u{0) = v'{0) = 1, u'{0) = v{0) = 0, whose 
Wronskian reads u{x)v'{x) — u'{x)v{x) = 1. We thus obtain 

(2.12) 



p2t,(L) + ip{u{L) + v'{L)) - u'{L) 
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The hypothesis of a smoothly varying potential implies that it is legitimate to use the 
well-known W.K.B. approximation PUJE]. Indeed, the condition for this scheme to be 
valid, 

K{xy dx ^ K{x)L ' ^ ' ' 

is automatically satisfied in the regime ()1.3|) of long barriers. Within this framework, a 
basis of solutions to (12.111 in the barrier reads 



i)±{x) ~ exp K{y) dy^ (0 < x < L). (2.14) 

At least one of the elementary solutions u{x) and v{x) (and generically both of them) is 
proportional to the growing solution ip^{x). Hence ()2.12p leads to the estimate 

(2.15) 



which will be sufficient hereafter, in order to work with exponential accuracy. 
The transmission therefore again falls off exponentially: 

InT ^ -2 K(x)dx, (2.16) 
Jo 

where the integral is the action of the classical imaginary-time trajectory, or instanton, 
crossing the barrier at energy E |1U1 ITTl IT^ . The estimate (|2.16|) generalizes (j2.1(jp to a 
potential barrier with an arbitrary (smooth) profile. 



3 Tunneling through a square disordered barrier 

We now turn to the case of a square disordered barrier of length L. The barrier potential, 

V{x)=Vo + W{x) (0<x<L), (3.1) 

is the sum of a constant Vq and of a disordered component W{x) with zero mean. 

As stated in the Introduction, in the regime ()1.3p it is legitimate to model W{x) as a 
Gaussian white noise, such that 

{W{x)W{y)) = 2DS{x -y). (3.2) 

The estimate ()2.15|) still holds in the presence of disorder. We are therefore led to 
consider the situation where the random potential extends over the half-line x > 0, and 
to investigate the growth rate of the generic solution to 1)2.11) . 

The most salient effect of the disordered potential is that the wavefunction in the 
barrier now changes sign many times in the regime under consideration, so that neither 
the concept of a single instanton trajectory, nor the W.K.B. estimate ()2.14|1 . make sense 
anymore. The most efficient approach to this problem is the invariant-measure method, 
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initiated long ago by Dyson and Schmidt [M]. For technical reasons, the energy E—iO 
and the inverse penetration length K + lO are respectively endowed with infinitesimal 
negative and positive imaginary parts. The method consists in introducing the Riccati 
variable 

= ^> (3-3) 

so that 

rx 

^jj{x) = ■ip{0) exp / z{y)dy. (3.4) 



The Schrodinger equation (|2.H) is equivalent to the Riccati equation 

z' = - + W{x). (3.5) 

The key property of this equation is the following: the complex random variable z{x) has 
a well-defined limit probability distribution, irrespective of the position x, provided it is 
deep enough in the sample {Kx ^ 1), and of the initial condition -2(0). Let us denote 
averages with respect to this invariant measure as ((...)). Equation (j3.4j) implies that the 
growing solution to (j2.ip typically grows exponentially, as 

ln?/;+(x) ^ Qx, (3.6) 

where 

Q = {{z)) (3.7) 

is the complex characteristic exponent. When the energy variable is just below the real 
axis, Q splits according to [121 E] 

n{E -lO) = -f{E) + 171 H{E). (3.8) 

The real part of ()3.8|) is the Lyapunov exponent 7. Inserting the behavior ()3.6|) into 
the estimate ()2.15|) . we therefore obtain at once the prediction 

(lnT)^-27L, (3.9) 

recalled in the Introduction [see ()1.2|) ]. The imaginary part of ()3.8|) is proportional to 
the integrated density of states of the problem per unit length, 

H{E) = r p{E')dE\ (3.10) 

J —00 

so that 1/H{E) is the mean distance between any two consecutive zeros of ip+i^)- 

The present case of a white-noise potential has been investigated by several au- 
thors jini ^1 ^] • Their findings can be recast into the following formula for the charac- 
teristic exponent: 

n = D^^^F{X), (3.11) 

with 
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and 

F(X) - e-^'W3 Ai-(e-^-/^X) _ Ai-(X) + i Bi-(X) 
^ " Ai(e-2W3X) " Ai(X) + iBi(X) ' ^^'^^^ 

where Ai(2;) and Bi{z) are Airy functions jil9j. As a consequence of ()3.8|) . we have 

7 = ^R(^) ^R(^) = Ai(X)Ai-(X) + Bi(X)Bi-(X) 
7 i7 ^ ^^j, ^ ^^j Ai(X)2 + Bi(X)2 ' ^""-'^^ 

ni/3 1 

References [13 HEl UZl HE] deal separately with the real and imaginary parts of Q, and 
therefore rather derive ()3.14|) and/or ()3.15|) . For completeness, we give in the Appendix 
a self-contained derivation of ()3.11|) - ()3.13|) . 

Inserting ()3.14|) into ()3.9p leads us to the prediction 

(InT) ^ -2D^/^ F^{X) L, (3.16) 

where the scaling variable X is real and positive in a tunneling situation. 

In order to emphasize the dependence of the transmission on the disorder strength, 
we recast the above prediction as 

(lnT)^(lnT)oG'(F), (3.17) 

where (InT)o = —2KL is the result ()2.10|) in the absence of disorder, while 

is the reduced disorder strength, and the scaling function G reads 

Weak-disorder regime 

The weak-disorder regime corresponds to D ^ K^, i.e., X ^ 1 or F ^ 1. The differential 
equation 

F^ + F' = X (3.20) 
obeyed by the function F{X) easily yields the asymptotic expansion 

F(X) = X^/2_J ^ + ---, (3.21) 

V J 4^ 32X5/2 ' V ; 

This result is formally real, so that F^ has the same asymptotic expansion (while F^ is 
exponentially small as X ^ +oo). We thus obtain 

Y 

G'(r) = l---^ + ---, (3.22) 

i.e., more explicitly 

{\nT)^2L(-K + ^ + ^^--]. (3.23) 

^ ' \ AK^ 32^5 J ^ > 

The first correction of order D is in agreement with the perturbative result of the 
leading effect of a weak disorder is found to enhance transmission. 
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Strong-disorder regime 

Consider now the opposite regime of strong disorder {D ^ K^, i.e., F 1 or X <^ 1). 
In this regime, the scaling function 

G{Y) ^ F^{0)Y^/^ (3.24) 

grows as the power 1/3 of the strength of disorder, with the exphcit prefactor 

^^(0) = ^ = 0.364506. (3.25) 

As a consequence, we have 

(InT) ^ -2F^(0) D^^^ L. (3.26) 

The transmission decreases with the disorder strength in the strong-disorder regime, as 
anticipated in the Introduction. The result ()3.26|) is independent of K: the damping due 
to the mean barrier Vq becomes negligible with respect to localization effects. 



Non-monotonic crossover behavior 

The function G(Y) which enters the scaling law ()3.17|) decreases first from its value 
Y{0) = 1, according to ()3.22|) . in the weak-disorder regime, and then increases according 
to ()3.24j) . in the strong-disorder regime. It must therefore reach a minimum, somewhere 
in the crossover from weak to strong disorder. 

Figure HI shows plots (full lines) of the scaling function G{Y), given by ()3.14|) . ()3.19|) . 
This function passes through a minimum, G{Y*) = 0.7284 for Y* = 1.695, before it 
crosses the value G{Y**) = 1 for Y** = 14.168. The non-monotonic behavior of the trans- 
mission as a function of the disorder strength in the crossover regime, anticipated in the 
Introduction, is therefore confirmed at a quantitative level. The tunneling transmission 
is enhanced for a weak enough reduced disorder strength {Y < Y^), the enhancement 
being maximal for Y = Y*. 

4 Tunneling through an arbitrary disordered barrier 

We finally turn to the general situation where both the deterministic part Vo{x) and the 
strength of the disordered part W{x) of the potential may have a smooth dependence on 
the position across the barrier. We set 

V{x) = Vo{x) + W{x) (0 < X < L), (4.1) 

with 

Kix) = ^Vo{x) - E, (4.2) 

and 

{W{x)W{y)) = 2D{x)6{x -y). (4.3) 
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Figure 1: Plots of the scaling functions describing the effect of disorder on tunneling transmis- 
sion, against the reduced disorder strength Y. Full lines: G entering l|3.17() (square barrier). 
Dashed lines: Gi entering (|4.9() (parabolic barrier with parabolic disorder). Dash-dotted lines: 
G2 entering (|4.14|) (parabolic barrier with uniform disorder), (a): Circles show the values of 
Y where the scaling functions cross unity: Y** = 14.168, Y^* = 11.469, Y2* = 6.522. (b) 
(enlargement): Circles show the values of Y where the scaling functions have their minima: 
G = 0.7284 for Y* = 1.695, Gi = 0.7331 for Y^ = 1.408, G2 = 0.7704 for Y^ = 0.982. 



We assume that K{x) and D{x) vary smoothly across the barrier, and we again 
focus our attention onto the regime ()1.H|1 . In this situation, the Riccati variable z{x) 
is approximately distributed according to the local invariant measure, characterized by 
the parameters K{x) and D{x). The reason why this adiabatic approach is legitimate 
is the same as for the W.K.B. scheme in the absence of disorder: the length over which 
parameters vary, i.e., the barrier length L itself, is much larger than the characteristic 
length of the relaxation of the distribution of the Riccati variable, i.e., the localization 
length 1/7(0;). 

Equations ()2.15p . ()3.4|) yield the following expression 

(InT) ^ -2 f^-f{x)dx (4.4) 
Jo 

for the mean logarithm of the transmission. The prediction ()4.4|) includes all the previous 
results fl2.10|) . ()2.16|) . ()3.9p as special cases. By means of ()3.14p . it can be recast into the 
more explicit form 

(luT)«-2^'z)(x)'/'f'(^)d.. (4.5) 

In the weak-disorder regime, this prediction behaves as 

(l„T>«2f (-A-W + ^ + ...)d.. (4.6) 
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while in the strong-disorder regime we have 

(InT) ^ -2F^(0) D{xfl^dx. (4.7) 
Jo 

These two expressions respectively generalize ()3.23p and ()3.26|) . 

To close up, let us illustrate more explicitly the result ()4.5|1 on two examples. 

Example 1: Parabolic barrier with parabolic disorder 

In our first example, both the deterministic potential and the disorder strength have a 
parabolic shape: 

K{xY = V,{x)-E = Kl^-^^^jpl, D(x) = Do^^4^, (4.8) 

with maximum values Kq and Dq at the center of the barrier (x = L/2). 
In this situation, ()4.5|1 yields the scaling law 



similar to (j3.17|) . with 



(lnT)^(lnr)oGi(F), (4.9) 



{\nT)^ = ~KoL, (4.10) 

y = ^.^ (4.11) 



and 



G^iY) = rF^fr-2/3(sin^)2/3) (sinOf^ AO. (4.12) 

The latter expression is obtained by setting x = L(l + cos6')/2, so that 4a;(L — x)/ L"^ = 
sin^ 9, with < 6* < vr. 

The scaling function Gi{Y) has been evaluated numerically by means of the inte- 
gral ()4.12|) . and plotted in Figure ^ (dashed lines). Its qualitative dependence on the 
disorder strength Y is similar to that of G{Y), with the following behavior at weak and 
strong disorder: 

Gi(r)^i-- iY<^i), 

1/3 ^^-^^^ 

G,{Y) ~ l^i5^lZ^yi/3 ^ 0.390454 yi/=^ {Y » 1), 

with a minimum Gi{Y*) = 0.7331 for Y^* = 1.408, and with Gi{Y-^**) = 1 for Y^** = 
11.469. 
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Example 2: Parabolic barrier with uniform disorder 



In our second example, the deterministic potential still has the parabolic form (j4.8|) . 
while the disorder strength D{x) = D is constant. 
The prediction ()4.5p now yields the scaling form 

(lnT)^(lnT)oG'2(r), (4.14) 

again with (|4.1U|) . and with 

Y = ^^ (4.15) 



and 



G^iY) = - F^/^ r ( sin' e) sin 6 dO. (4.16) 

The scaling function 6*2 (^) is also plotted in Figure [U (dash-dotted lines). Its quali- 
tative dependence on the disorder strength Y is again similar to that of G{Y), with the 
following behavior at weak and strong disorder: 

G2{Y) ^l-^\n^ iY<^l), 

(7) 

G,{Y) ^ 2 3^/^r(2/3) ^y, ^ Q 4g4^Q3 Yi/3 » 

71 r(i/3) 

with a minimum ^2(^2') = 0.7704 for Y* = 0.982, and with G2{Y^*) = 1 for Y** = 6.522. 



5 Discussion 



In this paper, we have investigated by analytical means the transmission through a disor- 
dered tunnel barrier, in the regime ()1.3|1 of most physical interest. We have thus provided 
a quantitative analysis of the phenomenon of disorder-induced enhanced tunneling, put 
forward by Freilikher et al. [HI IZ|- The most salient outcome of the present work is 
that the enhancement effect is a non-monotonic function of the strength of disorder, and 
that it is maximally efficient at some well-defined intermediate value Y* of the reduced 
disorder strength Y. 

The key point of our approach consists in utilizing the scaling law (j3.1ip - (|3.13p for the 
complex characteristic exponent Q. It is worth noticing that this formalism encompasses, 
and treats on the same footing, both the usual situation of the Anderson localization, 
where the energy is above the mean of the disordered potential (corresponding to neg- 
ative values of the scaling variable X), and the tunneling situation, where the energy 
is below the mean of the disordered potential (corresponding to positive values of X). 
As a consequence, the general results on the statistics of the transmission, recalled in 
the Introduction, still hold in the tunneling situation. In particular, (InT) is the right 
quantity to consider. 

In the case of a white-noise potential, considered throughout this work, the form ()3.12j) 
of the scaling variable X is merely dictated by dimensional analysis, while the explicit 
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formula ()3.13p for the scaling function F{X) is given a self-contained derivation in the 
Appendix. Somewhat equivalent results had been obtained in several earlier works 
HTHIH]. The scaling law (j3.1ip - (j3.13p also describes the spectra of other one- dimensional 
disordered systems, such as the diffusion of classical particles in a random force field |2D]. 
An analogous scaling formula holds for discrete models near their band edges. Consider 
the tight-binding equation 

Ipn-l + V'n+1 + Vn^pn = Slpn- (5.1) 

For a clean chain [Vn = 0), the dispersion relation reads S = 2 cos p. In the regime of a 
weak disorder {{V^) = A ^ 1), and near the upper band edge (p <^ 1), one has [21] 

This scaling law turns out to play a central role in various problems, such as the spreading 
dynamics of a wave packet [22]. The results ()3.1H) - ()3.13|) can be viewed as the formal 
continuum limit of ()5.2j) . The identification between our continuum problem and the 
tight-binding model (j5.ip has to take place near the band edge. Indeed, introducing ex- 
plicitly the lattice spacing a, which plays the role of the correlation range of the potential, 
the first inequality of ()1.3|) implies \pa\ = Ka <^ 1, as p = iK, whereas D = A/(2a^). 
The scaling variables and functions match between ()5.2|) and ()3.11|) - ()3.13|) (up to powers 
of 2 due to different conventions): x = 2^/^X, / = 2^/^F. 

For a square disordered barrier, the main result (j3.17|) for (InT), i.e., the logarithm of 
the typical transmission, can be generalized to all the moments of T in the regime ()1.3|) . 
It can indeed be shown, along the lines of |2] , that these moments scale as 

(T") ^ exp {-2n G„(F) KL) , (5.3) 

where Y is the reduced disorder strength of ()3.18p . while the scaling function G'„(F) 
depends on the order n (not necessarily an integer). Equation (|3.17|) is recovered as 
Go{Y) = G(Y). Skipping every detail, let us mention the weak-disorder expansion 

G„(r) = i ^Y + --- (r«i), (5.4) 



hence 



(T") ^ exp 




2L -nK+^-—^ + ... . (5.5) 



In agreement with \Z,, both (T) and (1/T), and indeed all the moments whose order is 
not in the range — 1/2 < n < 0, are increased by a weak amount of disorder. In the 
opposite strong-disorder regime, we have 

G„(F)^a„Fi/3 (>^>1), (5.6) 

with some n-dependent amplitude a„, hence 

(T") ^ exp (-2n a„ D^/^l) . (5.7) 



12 



The evaluation of the full scaling function Gn{Y) is a difficult task in general, except 
for n a negative integer, corresponding by means of (j2.15j) to positive integer powers of 
|'?/'+(L)p. In this case, Gn(Y) can be derived by means of the algebraic approach of |2j: 
it is (the real part of) an algebraic function with degree 2\n\ + 1. 

For a disordered barrier with an arbitrary profile, the general prediction ()4.5p has 
been illustrated on two realistic examples of parabolic barriers. These examples demon- 
strate that the qualitative features of disorder-induced enhanced tunnehng, and chiefiy 
its non-monotonic behavior as a function of the disorder strength, are rather insensitive 
to the shape of the mean and/or the width of the random potential. Even quantitative 
characteristics, such as Y* or Y**, do not depend too much on the profile of the barrier. 

Besides the present situation of thick tunnel barriers, the non-monotonic enhancement 
of transmission may be a more general phenomenon. Somewhat similar features have 
indeed been observed recently |23j in a problem inspired by nuclear fission. 

It is also worth noticing that, in the somewhat dual case of the total refiection of 
an electron by a semi-infinite disordered sample, universal features of the distribution of 
the Wigner time delay has received much attention recently j21] . The effects of disorder- 
induced enhanced tunneling on the Wigner time delay are also potentially of interest 
in transmission even though this concept has been criticized as being somewhat 
ambiguous ^Ej . 

Finally, the disorder-induced enhancement of transmission is a quantum-mechanical 
phenomenon, and more generally a wave phenomenon. It is indeed due to the existence of 
non-trivial localized states in a weakly-disordered barrier, which themselves originate in 
the interferences between the multiply scattered waves. Therefore, no disorder-induced 
enhancement is to be expected e.g. in the Kramers problem |2Z| of the thermally activated 
hopping of a classical particle over a potential barrier. 
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A Derivation of expressions ( 13.111) — (13.131) 



This Appendix presents a self-contained derivation of expressions (l3.1H) - (j3.13j) for the 
complex characteristic exponent fl, introduced in ()3.6p . 

To do so, we determine the invariant measure of the complex Riccati variable z, 
introduced in ()3.3p . As this distribution has a complex support, instead of writing a 
Fokker-Planck equation for its density, it is advantageous to consider the linear transforms 

= (iHy - z))) , = = • (A.l) 

For definiteness, we assume that K has positive real and imaginary parts. Equa- 
tion (|3.5|) shows that z{x) keeps a positive imaginary part, so that is analytic in the 
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lower half plane. 

In the spirit of the derivation of the Fokker-Planck equation |2H] , we consider a small 
increment e of x, and set z{x + e) = z{x) + rj. Both 

(r^) ^ (K^ - z^)e, {rf) ^ 2De (A.2) 

are proportional to the increment e, while higher cumulants are negligible. 

Along the lines of the Dyson-Schmidt approach jTSllEI, we then look for a stationarity 
condition by comparing the expressions of $(?/) corresponding to the points x and x + e: 
{{ln{y — z))) = ((ln(?/ — z — r]))). By expanding the right side of this equality in powers 
of 1], and using ()A.2|) . we obtain the condition 

i.e., 

D<p'{y) + {y' - K')<p{y) =y + Q. (A.4) 

The strength of disorder D can be scaled out by setting = D^/^X, Q = D^/^F, 
proving thus the scaling formulas ()3.11|) . ()3.12j) . and y = D^^^u, = D^^^^ip. 
We are thus left with a differential equation for iplu), 

iIj'{u) + {u'^-X)iIj{u)=u + F, (A.5) 

which can be solved by varying the constant: 

^Ij{u) =exp{-uy3 + Xu)C{u), C{u) = J {v + F) exp{vy3 ~ Xv) dv. (A.6) 

The existence of a regular solution, such that ip{u) — > as |m| ^ oo, determines F{X). 
We must have C{u) — as |m| — ^ oo in all the directions of the lower half plane where 
exp(— diverges. This happens in two Stokes sectors, represented by the directions 
u — > — oo and u — >• e~^'^/^oo. We thus obtain 

F - (A 7) 

with 

/g-i7r/3g(-j 
exp(t;V3 - Xv) dv. (A.8) 
-oo 

Finally, I{X) is equal to the Airy function Ai(e"^"^/^X) \19l, up to a multiplicative 
constant. This observation leads to ()3.13|) . 



14 



References 

[1] See e.g. I.M. Lifshitz, S.A. Grcdcskul, and L.A. Pastur, Introduction to the Theory of 
Disordered Systems (Wiley, New- York, 1988); B. Kramer and A. MacKinnon, Rep. 
Prog. Phys. 56, 1469 (1993); Y. Imry, Introduction to Mesoscopic Physics (Oxford 
University Press, Oxford, 1997); K. Efetov, Supersymmetry in Disorder and Chaos 
(Cambridge University Press, Cambridge, 1997). 

[2] J.B. Pendry, Adv. Phys. 43, 461 (1994). 

[3] Y. Imry and R. Landauer, Rev. Mod. Phys. 71, S306 (1999), and the references 
therein. 

[4] A. Crisanti, G. Paladin, and A. Vulpiani, Products of Random Matrices in Statistical 
Physics (Springer, Berlin, 1992). 

[5] J.M. Luck, Systemes desordonnes unidimensionnels (in French) (Collection Alea- 
Saclay, 1992). 

[6] V.D. Freilikher, B.A. Liansky, I.V. Yurkevich, A. A. Maradudin, and A.R. McGurn, 
Phys. Rev. E 51, 6301 (1995). 

[7] V.D. Freilikher, M. Pustilnik, and I.V. Yurkevich, Phys. Rev. B 53, 7413 (1996). 

[8] J. Heinrichs, Phys. Rev. B 33, 5261 (1986); 36, 2867 (1987). 

[9] L.I. Deych, A.A. Lisyansky, and B.L. Altshuler, Phys. Rev. Lett. 84, 2678 (2000). 

[10] See e.g. L.D. Landau and E.M. Lifshitz, Quantum Mechanics: Non-Relativistic The- 
ory (Pergamon, Oxford, 1959). 

[11] See e.g. N. Froman and P.O. Froman, Physical Problems Solved by the Phase-Integral 
Method (Cambridge University Press, Cambridge, 2002). 

[12] See e.g. J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Clarendon, 
Oxford, 1989). 

[13] F.J. Dyson, Phys. Rev. 92, 1331 (1953). 

[14] H. Schmidt, Phys. Rev. 105, 425 (1957). 

[15] Th.M. Nieuwenhuizen, Physica A 120, 468 (1983). 

[16] H.L. Frisch and S.P. Lloyd, Phys. Rev. 120, 1175 (1960). 

[17] B.I. Halperin, Phys. Rev. 139, A 104 (1965). 

[18] P.L. Sulem, Physica 70, 190 (1973). 



15 



[19] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions (Dover, New 
York, 1974). 

[20] J. P. Bouchaud, A. Comtet, A. Georges, and P. Le Doussal, Ann. Phys. 201, 285 
(1990). 

[21] B. Derrida and E. Gardner, J. Phys. (France) 45, 1283 (1984). 

[22] S. De Toro Arias and J.M. Luck, J. Phys. A 31, 7699 (1998). 

[23] B.G. Giraud, K. Amos, S. Karataghdis, and B.A. Robson (work in progress). 

[24] C. Tcxier and A. Comtet, Phys. Rev. Lett. 82, 4220 (1999), and the references 
therein. 

[25] C.J. Bolton-Hcaton, C.J. Lambert, V.I. Fal'ko, V. Progodin, and A.J. Epstein, Phys. 
Rev. B 60, 10569 (1999). 

[26] R. Landauer and T. Martin, Rev. Mod. Phys. 66, 217 (1994). 

[27] See e.g. V.I. Mel'nikov, Phys. Rep. 209, 1 (1991). 

[28] N.G. van Kampen, Stochastic Processes in Physics and Chemistry (North- Holland, 
Amsterdam, 1992). 



16 



